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Abstract 

Properties of one-dimensional (ID) arrays of low Ohmic tun- 
nel junctions (i.e. junctions with resistances comparable to, 
or less than, the quantum resistance R q = h/e 2 w 25.8 kO) 
have been studied experimentally and theoretically. Our 
experimental data demonstrate that - in agreement with 
previous results on single- and double-junction systems - 
Coulomb blockade effects survive even in the strong tun- 
neling regime and are still clearly visible for junction resis- 
tances as low as 1 kfi. We have developed a quasiclassical 
theory of electron transport in junction arrays in the strong 
tunneling regime. Good agreement between the predictions 
of this theory and the experimental data has been observed. 
We also show that, due to both heating effects and a rel- 
atively large correction to the linear relation between the 
half-width of the conductance dip around zero bias voltage, 
V1/2, and the measured electronic temperature, such arrays 
arc inferior to those conventionally used in the Coulomb 
Blockade Thermometry (CBT). Still, the desired correction 
to the half- width, AV]y 2 , can be determined rather easily 
and it is proportional to the magnitude of the conductance 
dip around zero bias voltage, AG. The constant of pro- 
portionality is a function of the ratio of the junction and 
quantum resistances, R/R q , and it is a pure strong tunnel- 
ing effect. 

1 Introduction 

The effect of Coulomb blockade in ID arrays of normal tun- 
nel junctions can be used for absolute thermometry @-(|]. 
The properties of such arrays have been extensively inves- 
tigated both experimentally and theoretically In all 
these works arrays of high Ohmic junctions with the junc- 
tion resistances, Rj, higher than or of the order of the quan- 



tum resistance R q — h/e 2 « 25.8 kf2 have been studied. In 
this limit a theoretical description of the Coulomb blockade 
is well developed || and it has been successfully applied Q- 
H to explain experimental findings in the high temperature 
regime, ksT > Eq = e 2 /2C, where C is the capacitance of 
a single junction. 

From the practical point of view, large (N 3> 1) ID arrays 
used as thermometers are advantageous over the smaller 
ones, because of higher accuracy of temperature measure- 
ments. At the same time, increasing the number of junc- 
tions in the array obviously yields an increase of its total 
resistance i? t ot °c N. Since in practice it is desirable to 
avoid very large values of i?totj it appears natural to relax 
the condition Rj 3> R q and consider ID arrays of relatively 
highly conducting tunnel junctions with Rj of the order of 
R q or smaller. A natural way to avoid this is the parallel 
connection of several ID arrays, and this is used extensively 
in Coulomb blockade thermometry (CBT). To avoid a large 
total number of junctions, a more straightforward solution 
would be to decrease the resistance of each individual junc- 
tion. On the other hand, heating effects turn out to be much 
more pronounced for highly conducting junction arrays [pj. 
Hence, from this point of view, it is better not to decrease 
Rj down to very low values. 

The above considerations motivated us to investigate the 
interplay between Coulomb blockade and strong tunneling 
effects in ID arrays of normal metallic tunnel junctions. Do 
Coulomb blockade effects survive if the junction resistance 
becomes smaller than the resistance quantum? Both theory 
and experiment Jl^]-[|l5| give a clear positive answer 
to this question. An adequate theoretical approach which 
enables one to study electron transport in the strong tun- 
neling regime is well established ||] . This so-called qua- 
siclassical Langevin equation technique allows to proceed 
analytically and remains accurate at not very low tempera- 
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Here cpj(t) = ^ J*dt'Vj(t') is the effective phase and Rs 
is the resistance of the electromagnetic environment. The 
shot noise of the j-th junction depends on tfj as 



R 1 ,C 1 




^—ei,cos(^)+6 J sin(^). 
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Here are Gaussian stochastic variables with the following 
pair correlators: 



Figure 1: A ID array of normal metal tunnel junctions. 



tures and/or voltages 
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for i? 3 ~ i?q, and max[/csT, eV] ~ e 2 /2C otherwise. The 
condition (Jlj) implies that in the strong tunneling regime 
i?j -C i? q , which is of a primary interest for us here, the 
technique Jl6], |^] covers practically all experimentally acces- 
sible values of temperature and bias voltage. In several pre- 
vious publications the Langevin equation technique 
was applied to analyze the Coulomb blockade and strong 
tunneling effects in single junctions and SET transistors, 
where the I — V curves were derived at arbitrary tunneling 
strength. The results of this theoretical analysis turned out 
to be in good agreement with available experimental data 
12l ~ 
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The structure of the paper is as follows: In Section 2 we 
develop a theoretical analysis of Coulomb blockade effects in 
ID tunnel junction arrays in the strong tunneling regime. 
In Section 3 we present experimental results obtained for 
the arrays with resistances in the range ^1 — 23 kS7 and 
compare these results with our theoretical predictions. Our 
main conclusions are summarized in Section 4. 



2 Theory 

In order to theoretically study the behavior of low resistance 
tunnel junction arrays we are going to use the technique of 
quasiclassical Langevin equation developed in Refs. |l(], M. 
In this section we generalize the approach presented in g 
to the case of ID arrays of tunnel junctions. 

We consider an array of iV normal metal tunnel junctions 
in series (Fig. 1). The system can be described by the 
following Langevin equations H-p|: 



2c + R, 2e 



,...,N; 



N 



v x =qRs + J2^r~ts R 



3 = 1 



2c 



(2) 



(5) 

where V stands for the principal value. There exist no cor- 
relations between the noise terms from different junctions: 
(6(2)i(*i)6(2)3> = for % ^ j. 

Averaging Eqs. (|^) over the noise realizations (we will de- 
note this average by angular brackets), we obtain the general 
expression for the current in the array: 



N 



V{V X ) 



(A) 



N 

E 

3=1 



"27 



(G) 



where i?s = Ylj=i Rj 1S the total resistance of the array. 
These two equations define the I — V curve of the chain. 

The problem, now, reduces to the evaluation of the noise 
averages (£,j). In order to do this, we first exclude the cur- 
rent q from the equations (^|) and get the following ones for 
the phases ipf. 
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Then we define the small deviations of the phase and the 
shot noise from their average values Stfij= (pj — (ipj) and 
<5£j = — respectively. They obey the following equa- 
tions 
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Applying the Fourier transformation and solving the corre- 
sponding equations, we find 
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where Zj(uS) is the total impedance seen by the j-th junction 
and the functions a,jk(io) describe the mutual influence of 
the junctions on each other. Equation ([)]) applies to any 
array of tunnel junctions of any dimensionality. In the case 
of ID arrays from Eq. (0) we find 



R , 



Rt, 



1 - iuR 3 C 3 R s + E- 

k 



(10) 



iiuRkCk 



Here we do not present the expressions for the functions 
ajk(i<j) because, as we will see below, the contribution of 
corresponding terms turns out to vanish. We find 
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where the response function Kj(t) is defined as 
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and Ajk(t) are defined analogously. 

It is important to emphasize that Eqs. (^|) and ( |i"l| ) 
are not the explicit solutions for the phase, but the inte- 
gral equations. The variables 5£j on the right hand side of 
these equations depend on 5ipj through the sin(t/3 :) /2) and 
cos(tfj/2) terms in the shot noise (||). These integral equa- 
tions can be solved by iteration. Here we restrict ourselves 
to the first iteration and put Sipj = in the right hand side 
ofEq. @. 

The next step is to evaluate the average values ) , which 
enter the expression for the current (||). We make the fol- 
lowing approximation: 
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Here Vj = ( ) is the average voltage on the j-th junction. 



Making use of Eq. (|11]) we get 
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Here we note that the terms containing the kernels Ajt~{t) 
do not contribute to the result because there exists no cor- 
relation between the noise on different junctions. Now the 



current is expressed as follows: 
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Let us first assume that all the junctions in the chain are 
identical, i.e. they have the same resistance R = Rj and 
the same capacitance C = Cj. Then we get 
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The current (0), then, can be found exactly 
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In the limit T — > the final result of the above expressions 
reduces to that of our previous analysis Q . 

The differential conductance is given by the following 
equation 
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Now let us put i?s = an d consider the high temperature 
limit u <C 1. 

Then, in the first order in u, we find 
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Figure 2: The normalized half- width of the conductance dip 
at zero bias voltage as a function of (inverse) temperature 
u = Ti/[2irkBTRC]. At higher temperatures the curve is 
almost linear in u [cf. Eq. (22) in the text]. 



This result exactly coincides with that found for the high 
Ohmic junctions (l], [|. 

Equation (pl| ) is basic for the Coulomb blockade ther- 
mometry. The half- width of the dip in the dl/dV is related 
to the temperature as = 

5A39Nk B T/e §. This rela- 
tion is proven to be very accurate for the high Ohmic arrays. 
To estimate its accuracy in the strong tunneling limit, we 
use the expression (|l9|), put i?g = and numerically solve 
the following equation 
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Here V1/2 is the normalized half- width, V1/2 = Vi/2/^1/2,0- 
The solution of this equation is plotted in Fig. 2. 

In the limit of high temperatures (small u) we find 
V1/2 = 1 + 0.704u - 0.24u 2 + • • ■ 



= 1 + 0.112- 



0.006- 



(23) 
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The zero bias conductance can be obtained from Eq. (|l^ 
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Here G% = 1/Rs, and 7 is the Euler's constant. Note 
that the result (^JJ) is equivalent to the one recently derived 
within the framework of a linear response approach based 



on the Kubo formula |L7). At high temperatures (small u) 
we find from (Q) 
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Now we will study the effect of the junction asymme- 
try on the properties of the thermometer. If the junctions 
are not identical, the equations become more complicated. 
Here, we will consider only the high temperature limit. In 
this limit the function G(t) decays at short times of order 
h/k B T, and one can replace the response function Kj(t) in 
the integral (|l4|) by its short time expansion, which starts 
from the linear-in-time term, Kj(t) — > Kj(Q)t. Then we get 
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Here Vj = VRj/R^ and Vj 
and ( |l5| ) we find 
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Assuming that deviations, <5i?j and <5Gj, of the junction 
parameters from the reference values R and G are small, we 
expand Eq. ( |27| ) in powers of and 6Cj up to the second 
order. Then we get 
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We observe that 5rj is zero if the deviations of all the 
resistances are equal to each other. We also note that the 
terms linear in SRj vanish in the sum Y]j Srj and, as a 
consequence, in Eq. (p^). Both these properties reflect 
the fact that the half-width depends only on temperature if 
the resistances of all the junctions are the same. Now the 
correction to the half-width of the conductance dip can be 
obtained peturbatively in SRj and 5Cj. In the first non- 
vanishing order we get 
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Figure 3: The temperature dependence of the zero bias 
voltage anomaly, R/ AR, for two arrays with per-junction 
resistances of 1.2 kil and 23 kil (solid triangles and open 
circuits, respectively). The solid lines are linear fits to the 
experimental data. 
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3 Experiment 

To obtain suitable data for comparison between the theoret- 
ical predictions (presented in the previous section) and the 
experiment, we fabricated high conductance Al/AlO^/Al 
tunnel junction arrays by electron beam lithography and 
two-angle shadow evaporation techniques. As a substrate, 
we used nitridized silicon wafers. The number of junctions 
in the array, N, was twenty, and each junction had an area of 
about 0.025 /im 2 . Different high conductance samples with 
per-junction asymptotic resistances of 1 — 2 kJl, together 
with two samples with lower conductances (in the interme- 
diate regime) with per-junction resistances equal to 20 kil 
and 23 kQ, were made and measured. In order to decrease 
heating effects in the high conductance arrays at higher bias 
voltages, the islands between the junctions in most of the 
samples (albeit not for those shown in Fig. 3) were made 
sufficiently large with cooling bars (see, e.g., M) attached 
to them. The measurements were carried out in the tem- 
perature range 1.5 K ~ T ~ 4.5 K, i.e. at temperatures of 
liquid helium. To measure the temperatures as accurately 
as possible, we fabricated CBT sensors on the same sample 
stage in the vicinity of the samples to be measured. 

As it was already discussed above, Coulomb blockade - 
although weakened - is not smeared out completely even 
in the strong tunneling limit. The zero bias conductance 
of the array is always lower than its asymptotic value at 



high voltages. In general, the zero bias conductance at high 
temperatures can be written as 
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In the limit of strong tunneling, Rj <C i? q , from Eq. (|2 



we find A = A 



— N-l 3C(3) fl q _ 



- jV-1 ga 



strong — N R — 0.0185 N R 



, while 



in the opposite weak tunneling limit the approach based on 



the Master equations, 0, yields A = A weak = ^ (^=1) . 
Here, as before, R = Rj is the (per-junction) resistance of 
the homogeneous array at large bias voltages. 

In the intermediate regime, which is appropriate for most 
of the arrays measured in the experiment, one can conjec- 
ture that A = A strong + A wea h. Previously, this conjecture 
was verified for the specific case of SET transistors (N = 2) 

Within the first order in u = TijlukBTRC the inverse 
resistance enhancement at zero bias voltage, R/AR, where 
AR = R(V = 0) — R, can be easily derived from Eq. 
as: 
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The expression above has two characteristic features: lin- 
earity in T and dependence of its slope on the capacitance 
of the junctions in the array. In addition, this equation pre- 
dicts an offset which depends only on the number of junc- 
tions in the array, and on the ratio of the quantum and 
per-junction resistances. Below, while making comparison 
between the measured data and the predictions of Eq. (|3^), 
we will take the the weak tunneling correction into account, 
i. e, A — A s i ron g -f- A wea k. 

The value R/AR, measured for two tunnel junction ar- 
rays with N — 20 and at different temperatures, is displayed 
in Fig. 3. The asymptotic resistances of the samples were 23 
kfi (open circles) and 1.2 kO (solid triangles). One observes 
an almost perfect linear dependence of the value R/AR on 
temperature (cf. Eq. (^3|). This dependence can be used to 
obtain a quantitative estimate for the junction capacitance. 
By fitting the slope of the experimental curves of Fig. 3 to 
Eq. © we find C = 2.4 fF and 2.1 fF respectively for the 
arrays with R = 23 and 1.2 Kf2. This way appears to be 
the most reliable to evaluate the per-junction capacitance 
in arrays of normal metal tunnel junctions Jl], 0]. 

The second and third terms in the right hand side of Eq. 
53]) yield an offset which - for a given number of junctions 
N - should depend solely on the ratio R q / R. Eq. ( |33| ) 
gives the values 3.4 and -0.2 for this offset respectively for 
the samples represented by the solid triangles and the open 
circles. The corresponding numbers obtained from the lines 
fitted to the experimental data in Fig. 3 are 4.1 and 0.2, 
being in a reasonable agreement with the above theoretical 
values. 

The offset values for two additional samples with (per- 
junction) resistances R = 2.1 kfi and R = 2.2 kfl (the sam- 
ples for which the systematic T-dependence was measured) 
were equal to 3.6 and 2.9, respectively. The corresponding 
theoretical predictions are 1.8 and 1.2. The data for the 
offsets are presented in Fig. 4. We observe that arrays with 
lower resistances show larger offsets, as predicted by our 
theory, Eq. (33). However, the offset values presented as 



a function of R q /R do not exactly fall on the straight line: 
they are scattered within an interval of ±0.5. This effect 
could probably be attributed to the inhomogeneity of the 
arrays. For instance, the measured offsets for two arrays 
with nearly identical values of R = 2.1 kS7 and 2.2 kSl dif- 
fer from each other. This discrepancy can be explained, if 
one assumes a certain degree of junction asymmetry, e.g., a 
20 % fluctuation of the per-junction resistance, Rj, around 
the mean value, R. This, in turn, can arise from a 40 % 
fluctuation in the area of the "odd" and the "even" tun- 
nel junctions as a consequence of the two-angle evaporation 
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Figure 4: The dependence of the offset of R/AR (cf. Fig. 
3) on the dimensionless resistance, R q /R, of the junctions. 
The measured data are shown by squares, the solid line is 
the best linear fit. 



technique employed in the sample fabrication. 

The best linear fit for the offset - inverse resistance de- 
pendence (Fig. 4) yields a = 0.19 ±0.05, b = 0.4 ±0.7. The 
value a agrees well with our theoretical prediction (|34|). At 
the same time the experimental value of b turns out to be 
different from the predictions of the theory (^). It is also 
too uncertain due to the scattering of the data points. At 
present, possible reasons for this discrepancy remain un- 
clear. We have checked asymmetry effects as well as those 
of the external environment within the framework of the 
Eq. (|l9|). These effects can hardly help to improve the 
agreement between theoretical and experimental values of 
b. 

In the weak tunneling regime R 3> R q the offset is given 
by b because aR q /R — > 0. In order to check the effect of 
finite external resistance on the offset in this limit we have 
performed Monte-Carlo simulations (for details see 0, |lS|| ) 
based on the phase-correlation theory |^0| . The only fitting- 
parameter used in these simulations was the junction ca- 
pacitance. By comparing the measured conductance curves 
of the array represented by the solid triangles in Fig. 3 
(R = 2.1 MY) to those obtained from the numerical simula- 
tions, we have obtained the value C = 2.1 fF for this sample. 
Notice that this value is in excellent agreement with that de- 
rived from Eq. ( ^3| ) as explained above. We have also found 
that, depending on the magnitude of the environmental re- 
sistance, b may increase by 0.1 at most. Again this is not 
sufficient to explain the observed values of b. 

Figure 5 shows the variation of the normalized half- 
width AV / i/ 2 /Vi/2 i o as a function of temperature (AV1/2 = 
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Table 1: Data for different samples at T ~ 4.2 K. Sam- 
ples with R = 2.1 kfi were further measured at lower tem- 
peratures (Fig. 5). The fifth row shows the measured 
data, whereas the last row has been obtained from the 
theory [Eqs. (36) and (37)]. Here we defined AVffi- = 



V, 



1 /2,corr 



^1/2,0 ( see a l so the text). 



Figure 5: The measured temperature dependence of the 
normalized half-width AVi/ 2 /V 1 / 2 .o for two arrays (solid 
and open squares) with the same per-junction resistances 
R = 2.1 kfi. The solid curve has been obtained from the 
strong tunneling theory (see the text). 

Vi/2— V\/2,o)- The solid and open squares represent the data 
measured for two 20-junction arrays both with R — 2.1 kfi. 
The solid curve corresponds to our theoretical prediction, 
Eq. (|23|). We observe, that the difference between the ex- 
perimental points and our theoretical curve typically does 
not exceed one percent, i.e. the ageement is fairly good. 

Within the framework of our analysis, one can easily es- 
tablish yet one more useful relation between the normalized 
half- width, A 1/1/2/^1/2,01 an d the value of the conductance 
dip AG = G s - G(V = 0) [cf, e.g., Eq. @] evaluated 
in the linear regime. Combining Eqs. ( f23| ) and (^5|) in the 
high temperature limit one gets 
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where the function x(Rq/ R) m the strong tunneling limit 



R q ^S> R reads 
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Note, that Eq. (|3^) was previously derived in Ref. 
the opposite weak tunneling limit R 3> R q , in which case 
the ^-function tends to the constant 



X{RJR) - 0.392. 



(38) 



We assume that in the intermediate regime the function 
x(Rq/R) is the sum of the expressions and (|38|). This 
assumption, although not proven rigorously, provides a rea- 
sonable interpolation between the two limiting cases. 



In addition to the arrays discussed above, three other 
samples with per-junction resistances of 20 kfi, 1.4 kfi, and 
1.0 kfi were measured at T ~ 4.2 K. The exact tempera- 
ture was obtained from the vapor pressure of liquid helium. 
The corresponding data are collected into Table 1. The sec- 
ond row in the table represents the measured half-widths 
obtained directly from the measurement, Vi/2,mcas- The 
corresponding conductance dips, AG/Gs, are shown in the 
third row. The forth row demonstrates the " corrected" half- 
widths (i.e. the weak tunneling correction [Eqs. ( |36| ) and 
(||) subtracted], 

y i/2 ,corr. = y/2,mcas.[l " 0.392^]. (39) 

In the weak tunneling regime this value would coincide with 
Vi/2,0 = 5.439A^fcsT/e. However in our experiments this is 
not the case. The fifth row shows the relative deviations 
of the corrected half-widths from that of the basic linear 
result, AVffio = (V 1/2 , C orr - V 1/2j0 )/V 1/2j0 . In fact, the 
numbers in this row can be interpreted as the "residual in- 
accuracies" of the measured half- widths as compared to the 
weak tunneling approximation. They can be explained by 
the strong tunneling effects. The corresponding theoretical 
values, obtained from the strong tunneling correction to the 
half-width of the conductance dip around zero bias voltage 
[Eq. ( [§6| , p7|)]) are shown in the last row of the table. These 
predictions are in a good agreement with the corrected mea- 
sured data discussed above. 



4 Conclusions 

We have studied one-dimensional arrays of tunnel junctions 
in the strong tunneling regime, both theoretically and ex- 
perimentally. Within the framework of the quasiclassical 
Langevin equation formalism, analytical expressions for the 
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current- voltage characteristics of such arrays, together with 
the expressions for the half-width of the conductance dip 
around zero bias voltage, have been derived. Furthermore, 
the effect of external elecromagnetic environment has been 
studied theoretically. We have fabricated and measured 
several arrays of tunnel junctions with per-junction resis- 
tances ranging from 1 kfi to 23 kf2. Our measured data are 
in rather good agreement with the theoretical predictions. 
The experiments demonstrate that Coulomb blockade ef- 
fects survive even in the strong tunneling regime, and are 
clearly visible in arrays with per-junction tunnel resistances 
as low as 1 kfi. These observations are in agreement with 
other recent experimental results on SET transistors |l2j 
and on single junctions |il|-||l5|. 

It has been verified, both in the theory and in the ex- 
periment, that high conductance arrays of tunnel junc- 
tion are less favorable for Coulomb Blockade Thermome- 
try (CBT) applications. This is a combined consequence 
of the more pronounced heating effects in such arrays, and 
the larger departure from the simple linear relation V1/2.0 = 
5.439A^fcsT/e, which is of central role in the CBT applica- 
tions. The latter, however, is not such a severe limitation, 
since corrections to the linear relation are now well known 
for all values of the junction resistance. It has also been 
shown that in the high temperature limit (in the leading 
approximation in l/T) the results for the conductance dip, 
now derived for arbitrary tunneling strength, coincide with 
those previously obtained in the weak tunneling regime. 
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